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in : Massive clusters of galaxies are the largest virialized objects in s Z e ffeVt;fo7a "review," see Sunyaev"& ZeY'dovich'l 980)," which 

O ■ the Universe bound by gravitation in the presence of dark mat- d dg Qn ^ details of ^ electron distribution function in the 

ter. Apart from galaxies and dark matter, galaxy clusters contain duster atmosphere (see> e Colafrancesco et al. 2003). 
a hot highly ionized plasma with temperatures up to 15 keV that 

emits in X-rays: a continuum through bremsstrahlung and lines A relativistically correct formalism for the SZ effect based 

through spontaneous decays of excitation states of ions. The hot on the probability distribution of the photon frequency shift 

X i intracluster plasma can be fully described by the particle distil- after scattering was given by Wright (1979) to describe the 

1 bution function (DF) Comptonization process of soft photons by mildly relativis- 

- C ?- 1 Equilibrium DFs are different in non-relativistic and rel- tic P lasma ' Relativistic corrections for the SZ effect allow 

ativistic statistics. A thermal equilibrium DF is described by us to measure the temperature of intracluster plasma (see, 

Maxwell-Boltzman and Juttner functions in the framework of e 'g" Pomtecouteau et al. 1998; Hansen et al. 2002) and have 

the non-relativistic and relativistic theories, respectively (Landau been studled both analytically (Colafrancesco et al. 2003, 

& Lifshitz 1976). High-energy phenomena allow us to study Colafrancesco & Marchegiam 2010) and by means of numeri- 

thermal mildly relativistic particle populations, such as mildly cal simulations (Prokhorov et al. 2010a). 

relativistic electron populations in clusters of galaxies. There are different methods proposed so far to derive a 

Inverse Compton (IC) scattering of hot free electrons in clus- plasma temperature by means of the SZ effect, which are based 

ters of galaxies on the cosmic microwave background (CMB) on the measurement of one of following quantities: the shift of 

radiation field is another effect that provides us with a method the crossover frequency (see Rephaeli 1995), the intensity slope 

to study hot plasmas in galaxy clusters, because IC scattering around the crossover frequency (see Colafrancesco et al. 2009), 

the wide frequency spectroscopy of the SZ effect spectrum espe- 

Send offprint requests to: D.A. Prokhorov e-mail: dally at high frequencies (Colafrancesco & Marchegiani 2010), 

phdmitry@stanford.edu and the ratio of the SZ intensities at two frequencies (Prokhorov 



Accepted . Received ; Draft printed: January 21, 2013 

ABSTRACT 

Aims. Measurements of the Sunyaev-Zeldovich (hereafter SZ) effect distortion of the cosmic microwave background provide methods 
to derive the gas pressure and temperature of galaxy clusters. Here we study the ability of SZ effect observations to derive the electron 
distribution function (DF) in massive galaxy clusters. 

Methods. Our calculations of the SZ effect include relativistic corrections considered within the framework of the Wright formalism 
and use a decomposition technique of electron DFs into Fourier series. Using multi-frequency measurements of the SZ effect, we 
find the solution of a linear system of equations that is used to derive the Fourier coefficients; we further analyze different frequency 
samples to decrease uncertainties in Fourier coefficient estimations. 

Results. We propose a method to derive DFs of electrons using SZ multi-frequency observations of massive galaxy clusters. We 
found that the best frequency sample to derive an electron DF includes high frequencies y=375, 600, 700, 857 GHz. We show 
that it is possible to distinguish a Juttner DF from a Maxwell-Bolzman DF as well as from a Juttner DF with the second electron 
population by means of SZ observations for the best frequency sample if the precision of SZ intensity measurements is less than 
0.1%. We demonstrate by means of 3D hydrodynamic numerical simulations of a hot merging galaxy cluster that the morphologies 
of SZ intensity maps are different for frequencies y=375, 600, 700, 857 GHz. We stress that measurements of SZ intensities at these 
frequencies are a promising tool for studying electron distribution functions in galaxy clusters. 

Key words, galaxies: cluster: intracluster medium; relativistic processes; cosmology: cosmic microwave background 

1. Introduction causes a change in the intensity of the CMB radiation toward 

clusters of galaxies (the Sunyaev-Zel'dovich effect, hereafter the 
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et al. 2010a). As was noticed by Rephaeli (1995) the correct rel- 
ativistic equilibrium distribution is essential for the proper inter- 
pretation of measurements of the SZ effect. 

One of the methods to study the relativistic equilibrium DF 
of particles is to perform fully relativistic molecular dynamics 
simulations (see Cubero et al. 2007 for the ID case; Montakhab 
et al. 2009 for the 2D case; and Peano et al. 2009 for the 3D 
case). However, to justify the use of the relativistic equilibrium 
DF in astrophysics, an observational confirmation is required. 

In this paper we propose a method to derive the relativistic 
equilibrium DF that is based on multi-frequency measurements 
of the SZ effect in massive clusters. This provides us with a 
more complete analysis of the electron distribution in the veloc- 
ity space than that given by knowing only the temperature value, 
and a method to verify the relativistic equilibrium DF. We also 
check the validity of using the Juttner DF as an appropriate ap- 
proximation to the universal (equilibrium) electron distribution 
in massive merging clusters. 

Hard X-ray emission tails reported in Beppo-SAX and 
RXTE X-ray spectra of some galaxy clusters (see Fusco- 
Femiano et al. 1999; Fusco-Femiano et al. 2004; Rossetti & 
Molendi 2004, for the Coma cluster; Petrosian et al. 2006 for the 
Bullet cluster) were interpreted as bremsstrahlung emission from 
non-thermal subrelativistic electrons (e.g. Sarazin & Kempner 
2000) or from thermal electrons with a Maxwellian spectrum 
distorted by a particle acceleration mechanism (Blasi 2000; 
Liang et al. 2002), or from thermal electrons with a Maxwellian 
spectrum with a high temperature (e.g., Petrosian et al. 2006). 
Million & Allen (2009) reported the discovery of spatially ex- 
tended, non-thermal-like emission components in Chandra X- 
ray spectra for five massive galaxy clusters. Using Swift/BAT 
data, Ajello et al. (2010) have confirmed the presence of a hard 
X-ray excess from the Bullet cluster. Note that a possible con- 
tribution from inverse Compton emission of highly relativistic 
electrons with Lorentz factor y ~ 10 4 to a hard X-ray excess (for 
a review, see Rephaeli et al. 2008) constrains our possibilities 
to study electrons DFs by using X-ray observations. It has been 
suggested that a good test to check the bremsstrahlung interpre- 
tation of hard X-ray tails is to use multi-frequency measurements 
of the SZ effect (see Dogiel et al. 2007). 

Spatially-resolved observations of the SZ effect can also 
provide relevant information on the distribution of the electron 
plasma in galaxy clusters with various temperatures. Prokhorov 
et al. (2010a) found that the morphologies of the SZ intensity 
maps of a cool galaxy cluster at frequencies of 128 GHz and 369 
GHz are similar. Here we perform numerical simulations of a 
hot merging galaxy cluster to study the morphologies of the SZ 
intensity maps for a hot merging galaxy cluster to explore the 
different morphologies for low and high frequencies that depend 
on the importance of the SZ effect relativistic corrections. 

The layout of the paper is as follows. The dependence of the 
CMB distortion caused by the SZ effect on the electron DF is 
considered in Sect. 2 in the framework of the relativistic correct 
Wright formalism. We propose a method to derive a velocity DF 
from multi-frequency SZ observations in Sect. 3. Relaxation of 
a system of electrons to equilibrium distributions is considered 
in Sect. 4. We estimate the precision of SZ observations, which 
is necessary to derive the electron DF proposed as an explana- 
tion of hard tails in X-ray spectra of galaxy clusters, in Sect. 5. 
Using 3D hydrodynamic numerical simulations of a hot merging 
galaxy cluster and the Wright formalism, we demonstrate that 
the morphologies of the SZ intensity maps are different at low 
and high frequencies in Sect. 6. We present the observational es- 



timates and discuss the confusion noises in Sect. 7. We discuss 
our results and present our conclusions in Sect. 8. 

2. Dependence of the CMB distortion caused by the 
SZ effect on the DF of electrons 

In this section we discuss the difference between the CMB dis- 
tortions (caused by the SZ effect) that are caused by a departure 
from the diffusive approximation given by Kompaneets (1957) 
from those that are caused by using a relativistic correct DF in- 
stead of a Maxwell-Boltzman DF. 

In the diffusion approximation and for a non-relativistic elec- 
tron population, the CMB intensity change caused by the SZ ef- 
fect is 

kbT 

A/ nr = I q t y m (x) (1) 

where Iq = 2(£/,7cmb) 3 I (he) 2 , T e is the electron temperature, r 
is the the scattering optical depth, m e the electron mass, c the 
speed of light, kb the Boltzmann constant, x = hv/(kbTcMs), h 
the Planck constant, and the spectral function g m (x) is given by 

jcV / e x + 1 \ 

^ = W^W'—i " 4 I (2) 

The subscript 'nr' denotes that the previous expression was ob- 
tained in the non-relativistic limit. 

In the relativistic treatment taking into account the scattering 
to arbitrary frequencies (Wright 1979), the CMB spectral dis- 
tortion caused by the SZ effect is a functional of the electron 
velocity DF and is given by 

M[x, M = /or f dsPl [,,/.] - T^l) (3) 

with 

Pi[s,f e ] = J dpf e (0)P(s,0) (4) 

where P\ [s, / e ] is the probability distribution of a scattering fre- 
quency shift that is a functional of the electron velocity DF, f e (J3), 
of the electron population and P(s,/3)ds is the probability that a 
single scattering of a CMB photon off an electron with speed /3c 
causes a logarithmic frequency shift s = ln(v' /v), and p — V/c, 
where V is the electron velocity. 

We notice that this formalism is valid in the single scattering 
approximation and for low values of the optical depth r, which is 
however sufficient for the purposes of this paper. A more general 
description of the CMB spectral distortions caused by the SZ 
effect can be found in Colafrancesco et al. (2003). 

The SZ intensity spectra G(x) = AI(x)/(Ioy), where y is 
the Comptonization parameter, derived in the framework of the 
Wright formalism assuming Juttner and Maxwell-Boltzman DFs 
are shown in FigQ]for a massive cluster with a temperature of 
15.3 keV. The SZ intensity spectrum in the Kompaneets approx- 
imation is also shown in FigQ] 

To study the effects of the CMB intensity change caused 
by the SZ effect that is in turn caused by a departure from 
the diffusive approximation and by using a relativistic correct 
DF instead of a Maxwell-Boltzman DF, we calculate the ratios 
Al [fj(J3)] /A/nr and A/ [fu(fi)] /A/ nr as functions of the electron 
temperature, where fj(J3) and fyi(fi) are the Juttner and Maxwell- 
Boltzman DFs, respectively. These ratios are shown in Figj2]for 
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Fig. 1. SZ intensity spectra G(x) = AI(x)/(Ioy) for a massive 
cluster with a temperature of 15.3 keV for Juttner and Maxwell- 
Boltzman DFs shown by the solid and dashed lines, respectively. 
The SZ intensity spectrum in the Kompaneets approximation is 
shown by the dotted line. 

Juttner and Maxwell-Boltzman DFs by solid and dashed lines, 
respectively, at a frequency of 369 GHz, the choice of this fre- 
quency is reasonable since the relativistic corrections of the SZ 
effect are larger at high frequencies. 




T (keV) 

Fig. 2. Ratios A/ [/}(£)] /AI m and AI[f M (fi)] /A/ nr as functions 
of the electron temperature are shown by the solid and dashed 
lines, respectively 

To qualify the CMB intensity change caused by the SZ effect 
that is in turn caused by using a relativistic correct DF instead 
of a Maxwell-Boltzman DF, we calculate the ratio of the CMB 
intensity change caused by the SZ effect that is in turn caused by 
using a relativistic correct DF instead of a Maxwell-Boltzman 
DF to that given by the total contribution of the SZ relativistic 
corrections. This ratio is given by the expression 



D = 



M[f M m-M[fm] 



(5) 



and is shown (in %) at a frequency of 369 GHz in Figj3] as a 
function of temperature. Figure |3]shows that the correction from 
using a relativistic correct DF instead of a Maxwell-Boltzman 
DF is a more significant fraction of the total contribution of the 
SZ relativistic corrections at lower temperatures. However, we 
stress that the value of the CMB intensity change because of 
using a relativistic correct DF instead of a Maxwell-Boltzman 
DF will be much higher in hot galaxy clusters and, therefore, 



Fig. 3. Ratio of the CMB intensity change caused by the SZ 
effect caused by using a relativistic correct DF instead of a 
Maxwell-Boltzman DF to that given by the total contribution of 
the SZ relativistic corrections. 



measurements of the relativistic SZ corrections in hot clusters 
will be more effective to probe the electron DF. This is because 
the value of the relativistic SZ corrections is proportional to T^ 2 , 
since t oc sfTg (e.g., Bryan & Norman 1998) and 



J dsPi[s,f e ][ 



-3 s 



exp(x<? s ) - 1 e x - 1 



?n e c 2 



s nr (x)<xr e 2 ,(6) 



see, e.g., Challinor & Lasenby (1998). Therefore, the relativistic 
SZ corrections for a galaxy cluster with a gas temperature of 15 
keV will be in a; 16 times larger than that for a galaxy cluster 
with a gas temperature of 5 keV. Because the second-order and 
third-order (in the expansion parameter = ki,T e /(m e c)) rela- 
tivistic effects make a significant contribution to the SZ spectral 
distortion for k h T e =s 15 keV (Challinor & Lasenby 1998; Itoh 
et al. 1998), we checked the relativistic SZ corrections for a hot 
galaxy cluster by using the Wright formalism and found that the 
relativistic SZ corrections at frequencies of 369 GHz and 857 
GHz for a galaxy cluster with a gas temperature of 15 keV are 
«14 and «17 times larger than those at frequencies of 369 GHz 
and 857 GHz for a galaxy cluster with a gas temperature of 5 
keV, respectively. This agrees with the result obtained above. 

Because the the SZ intensity change from using a relativistic 
correct DF instead of a Maxwell-Boltzman DF is a significant 
fraction of the total contribution of the SZ relativistic corrections 
(see Figj3]l, we conclude that this should provide us with the 
ability to derive an electron DF by means of the SZ effect. In 
the following section we propose a method to derive a DF of 
electrons from multi-frequency SZ observations. 



3. A method to derive a velocity DF from 
multi-frequency SZ observations 

We use a Fourier analysis to derive a DF of electrons from multi- 
frequency SZ observations and to find the best frequencies at 
which this method can be successfully applied. 

We decompose the Juttner and Maxwell-Boltzman functions 
(in the beta representation /3 = V/c) into the Fourier cosine series 
and note that Maxwell-Boltzman and Juttner functions are ap- 
proximated with a high precision by six terms of this decomposi- 
tion (these terms are 1/2, cos(7r/?), cos(2tt/3), cos(3n/3), cos(4jt/3), 
and cos(57Ty6)). For a plasma with a temperature of 15.3 keV 
(which corresponds to = kbT e /(m s c) = 0.03), the decomposi- 
tion coefficients of a Juttner function approximately equal 2.000, 
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1.256, -0.172, -0.954, -0.786, -0.310 for the terms 1/2, cos(tt/3), 
cos(2jt/3), cos(3tt/3), cos(47ryS), and cos(57Ty6), respectively. For a 
plasma with a temperature of 15.3 keV, the decomposition co- 
efficients of a Maxwell-Boltzman function approximately equal 
2.000, 1.214, -0.204, -0.878, -0.699, -0.316 for the terms 1/2, 
cos(7ryS), cos(27t/?), cos(37Ty6), cos(47r/3), and cos(57ryS), respec- 
tively. Note that the coefficient at the first term equals 2, be- 
cause DFs are normalized. Juttner and Maxwell-Boltzman func- 
tions for a temperature 15.3 keV are shown in Fig J4] by solid 
and dashed lines, respectively. The approximation to a Juttner 
distribution obtained by means of six terms of Fourier series is 
shown in Fig. |4]by a dotted line. Figure |4] shows that six terms 
of the Fourier series is sufficient to describe the main features of 
DFs: an increase of the Juttner DF in the range /3=[0.2, 0.4] and 
a decrease of the DF in the range /3-[QA, 0.6] with respect to the 
Maxwell-Boltzman DF, which arise from relativistic corrections. 



fork=0, 1,2, 3,4, 5; and 



K(s, x) = 



exp(xe s ) - 1 e x - 1 



(10) 



Multiplying the left-hand and right-hand sides of Eq.© by 
l/0(r e ), we rewrite Eq.® in terms of the generalized spectral 
function G(x, T e ) in the form 



Al 



/oT®(r e ) 



= A k G(x, - 1 /26 k + cos (nkp)). 



(11) 



Since Aq = 2 because of the DF normalization, we need to 
find only five Fourier coefficients. Using multi-frequency mea- 
surements of the SZ effect, we can find a solution for the linear 
system of equations 




0.4 0.6 
beta 



A/fa) 

/oT®(r e ) 



J] A k Gfa, -l/26 k + cos (nk/3)) 



(12) 



k=0 



in order to derive the Fourier coefficients A k , for A=l, 2, 3, 4, 5. 

Note that the Comptonization parameter should be indepen- 
dently derived from a precise SZ intensity measurement at a fre- 
quency of v = 255 GHz (x=4.5), since the SZ intensities for 
Juttner and Maxwell-Boltzman DFs have the same value at this 
frequency (see FigJTJ. 

To make the system of linear equations well-conditioned, we 
constrain our analysis of DFs by considering the functions for 
which f e (J3 — 0) = 0. This choice is consistent with both Juttner 
and Maxwell-Boltzman DFs. In this case, taking into account 
that Aq = 2, the system of linear equations (fT2l can be written as 



Fig. 4. Relativistic and non-relativistic DFs, and the approxima- 
tion to a Juttner distribution obtained by means of six terms of 
Fourier series shown by the solid, dashed, and dotted lines, re- 
spectively. 



To show how a velocity DF of electrons can be derived from 
multi-frequency SZ observations, we write Eq.(|3} by using the 
generalized spectral function G{x, T e ) introduced by Prokhorov 
et al. (2010b) given by the expression 



0(r e ) \exp( xe~ s ) - 1 e x - 1 



ds. 



(7) 



Thus, the CMB intensity change caused by the SZ effect is 
Al = Iot (kbT e /(m e c 2 )} x G{x,T e ). Note that in this notation, 
the spectral function g m {x) is changed to the generalized spec- 
tral function G(x, T e ), which explicitly depends on the electron 
temperature. 

The CMB intensity change caused by the SZ effect for a DF 
represented by six terms of Fourier cosine series is 



Al = I t ^ A k f dsP x 0,-1 /26 k + cos (nkfi)] K(s, x), 

i,_n *J 



(8) 



where 6 k is the Kronecker delta function, Ak are the Fourier co- 
efficients given by 



Ai 



Jo 



f e (j3) cos(nkx)dx 



(9) 



J] A k (G(x A ,-l/26 k + cos (nk0)) - g m -fa)) = 



(13) 



rfa) + 



A/fa) 

/ r©(r e ) 



®(T e ) 



K(s, x A )ds. 



To check if this system of linear equations is well- 
conditioned, we study the properties of the matrix M Ak , which 
is given by 



Mm = Gfa, -l/26 k + cos (nk/3)) - g m (x A ) 



(14) 



for different samples of frequencies. 

To find a suitable sample of frequencies for deriving the elec- 
tron DF, let us consider the following samples: the first sample is 
in the low-frequency range including four frequencies v = 100, 
120, 140, 160 GHz considered by Colafrancesco & Marchegiani 
(2010); the second sample including four frequencies in the 
frequency range 300-400 GHz considered by Colafrancesco & 
Marchegiani (2010), v = 300, 320, 340, 360 GHz; the third sam- 
ple including low and high frequencies v = 100, 200, 300, 400 
GHz; and the fourth sample including more higher frequencies 
v = 375, 600, 700, 857 GHz. 

The first detection of the Sunyaev-Zel'dovich effect increment at 
such high frequencies where the relativistic corrections of the SZ 
effect are relevant has been obtained with HERSCHEL-SPIRE 
for the Bullet cluster (Zemcov et al. 2010). 

To find the frequency sample for which the electron DF can 
be optimally derived, we calculate the condition numbers for the 
matrixes of Mm and the inverse matrixes of MrJ. Since these 
samples have only four frequencies, the equation obtained from 
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the condition of f e (J3 = 0) = is used A\ +A 2 +A-}+A4+A 5 = -1 
(the normalization condition Ao = 2 is taken into account) to 
form the matrix 5 x 5, we choose the first row of these matrixes 
equaled to M\k = (1,1,1,1, 1). The condition numbers based on 
the L 2 norm of these matrixes approximately equal to 4.4 x 10 4 , 
2.1 x 10 5 , 1.6 x 10 3 , and 323.8 for the 1st, 2nd, 3rd, and 4th 
samples, respectively. Since the condition number of the matrix 
corresponding to the fourth frequency sample is the smallest, we 
conclude that the best sample to determine DFs by means of the 
SZ effect is the fourth sample, which includes high frequencies. 
The worst sample is the second sample, which is in the frequency 
range 300-400 GHz, since the matrix M^p) is ill-conditioned. 

To clarify our method of the DF analysis, we calculate the 
inverse matrixes 



where x\,X2,x-$,Xn correspond to the considered frequencies (v 
= 375, 600, 700, 857 GHz for the fourth sample) and the gen- 
eralized spectral functions G(x A , T e ) are calculated for a Juttner 
DF (this is because we study the ability of multi-frequency SZ 
observations to distinguish a Juttner DF from other DFs, which 
are slightly different from a Juttner DF). 
Equation d20b can be written as 



= 2 W lA V A 



where 



oA ■ 



(23) 



(24) 



M, 



-0.86 
-3.11 
-5.28 
2.00 
8.25 



0.15 
-19.35 
190.19 
159.65 
-330.66 



-0.59 
18.60 
-522.48 
-358.39 
862.87 



1.59 
5.30 
499.98 
294.83 
-801.71 



-1.05 
-6.74 
-164.28 
-89.09 
261.18 



(15) 



The matrixes W s ,\ are shown below for the samples of fre- 
quencies v = 375, 600, 700, 857 GHz and v = 300, 320, 340, 
360 GHz for a comparison. 



M, 



kA(2) 



-0.89 -6.89 
-3.83 42.01 
-6.68 808.77 
2.79 331.94 
9.63 -1175.84 



25.02 
-112.45 
-2462.02 
-1046.14 
3595.60 



-28.79 
97.69 
2509.77 
1104.08 



10.69 
-27.23 
-856.15 
-390.21 



-3682.75 1262.90 



Mi 



ki(3) 



M, 



M(4) 



-0.88 0.28 -0.14 -0.03 0.02 
-3.75 -0.90 1.76 -2.14 1.01 
-6.76 -22.04 25.24 -11.68 1.06 
2.57 -9.60 -1.86 -1.86 -1.71 
9.82 32.25 -36.24 15.70 -0.38 



-0.21 0.31 -2.18 2.79 -0.96 
-2.36 -0.06 -2.43 3.95 -1.61 
-2.79 -0.38 -2.06 5.24 -3.19 
0.86 -1.00 6.94 -9.18 2.95 
5.49 1.12 -0.28 -2.79 2.81 



(17) 



(18) 



The experimental uncertainties in the SZ intensity measure- 
ments constrain the ability to derive DFs of electrons. The un- 
certainties in the determination of the Fourier coefficients of At 
are given by 



SA J = Ta M ^ 

a 



6(AI(xd/h) 



and can be written as 

11 A 

where 



V., = 



( 6(M(x)) 
{ A7(x) 



(19) 



(20) 



(21) 



for A — 2, 3, 4, 5 and V\ — (since the first term determined from 
the condition f(J3 = 0) = 0), and 



G n A 



fQ 

G(x u T e ) 
G(x 2 ,T e ) 
6 G(x 3 ,r e ) 
























G(x4, T e ) 



(22) 



jyt(4) - 



(16) 



jA(2) 



I 





I 







1.75 
-0.33 
-2.10 
-5.54 

6.22 



-6.27 
-6.99 
-5.92 
19.99 
-0.80 



4.61 
6.53 
8.67 
-15.19 
-4.62 



-27.07 
165.02 
3176.87 
1303.88 
-4618.70 



115.26 
-517.97 
-11340.09 
4818.54 
16561.34 



-0.62 ^ 

-1.03 

-2.04 

1.89 

1.80 



-146.87 
498.23 
12799.83 
5630.83 
-18782.03 



(25) 



57.94 
-147.59 
-4640.34 
-2114.94 
6844.92 



(26) 



Since the absolute values of the elements of W-^a) are much 
smaller than those of Wj.jp), we justify our previous conclusion 
that the fourth sample is more suitable for the SZ analysis of 
electron DF than the second sample, since in the case of the 
fourth sample to derive a DF of electrons the allowable uncer- 
tainties in the SZ intensities are about three orders of magnitude 
larger (see Eqj23j. We also checked that the fourth sample is the 
best amongst the all samples considered in this paper in order to 
study DFs. 

To study the ability to derive DFs from SZ intensity mea- 
surements, we calculate the differences between the correspond- 
ing Fourier coefficients for Juttner and Maxwell-Boltzman DFs, 
which are \A U - A IM b\ = 0.042, |A 2 ,j - A 2 , M bI = 0.032, 
|A 3 ,j-A 3 ,mbI = 0.076,'|A 4 ,j-A 4 ,mbI = 0.0~87, and |A 5J -A 5 ,mbI = 
0.006 (see the previous section). 

Note that the difference \A^j - As mbI = 0.006 is very small 
and it will be difficult to observe this difference by the SZ mea- 
surements. However, we checked that this difference has a neg- 
ligible impact on the precision of the approximation of a Juttner 
DF obtained by means of Fourier series because of its small 
value. 

For the sake of illustration, we assume that the relative SZ 
intensity uncertainty in measurements d(AI(x,\))/AI Xl does not 
depend on the frequency and equals 

From Eqs.(l23l and d25l l we find that the uncertainty in SZ 
intensity measurements to distinguish the values of the Fourier 
coefficients for Juttner and Maxwell-Boltzman DFs should be 
less than £ = 0.32% to derive the A\ value, £ = 0.22% to derive 
the A 2 value, £ = 0.41% to derive the A 3 value, and £, = 0.20% 
to derive the A4 value. 

To check the result that an electron DF can be derived by 
means of the SZ effect if uncertainties of SZ observational data 
are less than 0.2 %, we performed Monte-Carlo simulations of 
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Fig. 5. Absolute difference between the Juttner and approxi- 
mate DFs for a temperature of 15.3 keV (solid lines) compared 
with the absolute difference between the Juttner and Maxwell- 
Boltzman DFs (dashed line). 



SZ observations (calculated in the relativistic correct formalism 
and using a Juttner DF) with the SZ intensity uncertainty of 0. 1 % 
at the frequencies included in the fourth sample. Using the SZ 
spectra obtained from the Monte-Carlo simulations, we solve 
the system of equations (13) to find the Fourier coefficients Ay_ 
and the approximate functions to a Juttner DF, which has been 
initially used to produce the CMB intensity changes caused by 
the SZ effect in Monte-Carlo simulations. The absolute value of 
the difference between the Juttner and Maxwell-Boltzman DFs 
for a temperature of 15.3 keV in the range of /3 — [0.2,0.6] is 
shown in Fig|5]by a dashed line. This range of /? contains the 
main features of these DFs: an increase of the Juttner DF in the 
range (3 = [0.2, 0.4] and a decrease of the Juttner DF in the range 
P = [0.4, 0.6] with respect to the Maxwell-Boltzman DF. The 
absolute difference between the Juttner DF for a temperature of 
15.3 keV and the approximate functions to the Juttner DF ob- 
tained by the Monte-Carlo simulations are shown in FigJ5]by 
solid lines. Figure [5] shows that the approximate DFs are close 
to a Juttner DF and it is possible to distinguish the Juttner and 
Maxwell-Boltzman DFs, if the SZ intensities are measured with 
an uncertainty of 0.1% (see Sect. 6 below). 

We conclude that the SZ effect provides us with an interest- 
ing method to study electron DFs in massive galaxy clusters that 
contain hot plasmas with temperatures (^15 keV) . 



4. Relaxation of a system of electrons to 
equilibrium distributions 

To justify the use of a Juttner distribution function as an appro- 
priate approximation to the universal (equilibrium) electron dis- 
tribution in hot merging clusters, we calculate the electron equi- 
libration time and compare this time with the merging time scale. 
We consider here the relaxation to equilibrium distributions in 
the framework of both the non-relativistic and relativistic theo- 
ries. 

We solve numerically the time-dependent Fokker-Planck 
equation for the evolution of an isotropic system of electrons 
with Coulomb interactions and with an initially Gaussian mo- 
mentum distribution. The Fokker-Planck equation coefficients 
are taken from Liang et al. (2002) and Dogiel et al. (2007). The 
initial Gaussian momentum distributions are chosen so that the 
number density normalization and the total energy are consis- 
tent with those of the equilibrium distributions. Note that both 



the Fokker-Planck coefficients (see Dogiel et al. 2007) are pro- 
portional to m e c 2 l(k\,T e ) and, therefore, the value of the char- 
acteristic frequency (which determines the rate of relaxation) is 
given by the inverse Spitzer collision time (see Eq. 18 of Liang 
et al. 2002). We use here the dimensionless time t in units of the 
Spitzer collision time. For an electron population with number 
density of n = 10~ 3 cm -3 and temperature T e = 15.3 keV, the 
Spitzer collision time equals 1.3 X 10 5 yr. 

The results of our numerical analysis of relaxation of a sys- 
tem of electrons to equilibrium distributions are plotted in Fig. 
i 




Fig. 6. Difference between the electron distribution function and 
the Maxwell-Boltzman DF at time t = 40 (solid) and 60 
(dashed) is shown as a function of {S = v/c. The analogous dif- 
ference between the electron distribution function and the Juttner 
DF at time t — 40 (dash-dotted) and 60 (dotted) is also shown. 



The solid and dashed curves in Fig|6]show the difference be- 
tween the electron distribution function, derived in the frame- 
work of the non-relativistic theory, at times r = 40 and 60, re- 
spectively, and the reference Maxwell-Boltzman DF. Since these 
differences are significantly smaller than the difference between 
the Juttner and Maxwell-Boltzman DFs shown in Fig. |5j the 
electron DFs at time r = 40 and 60 can be considered as good 
approximations to the Maxwell-Boltzman DF. 

The dash-dotted and dotted curves in Fig|6] show the dif- 
ference between the electron distribution function, derived in 
the framework of the relativistic theory, at times r = 40 and 
60, respectively, and the reference Juttner DF. These differences 
are again significantly smaller than the difference between the 
Juttner and Maxwell-Boltzman DFs shown in Fig. [5] and, there- 
fore, the electron DFs at a time of r = 40 and 60 can be consid- 
ered as good approximation to the Juttner DF. 

We stress that the qualitative estimate of the time required to 
fill the Maxwell tail to the velocity of 0.7c in the non-relativistic 
theory (see, e.g., MacDonald & Rosenbluth 1957) is twenty five 
Spitzer collision times and is in rough agreement with our nu- 
merical analysis. 

Since sixty Spitzer collision times correspond to « 8.1 x 10 6 
yrs for the number density and temperature values mentioned 
above, this time scale is a tiny fraction of the Hubble time. The 
merging time scale of a merging cluster equaled to 10 8 - 10 9 
yrs (see, e.g., Prokhorov & Durret 2007; Akahori & Yoshikawa 
2010) is longer than sixty Spitzer collision times. Thus, we con- 
clude that the Juttner distribution function is an appropriate ap- 
proximation to the universal (equilibrium) electron distribution 
in hot merged clusters. 
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We estimate the gas density at which the time-dependent DF 
relaxing to a Maxwell-Boltzman DF is indistinguishable from a 
Juttner DF. The time, at which the difference between the time- 
dependent DF and the Maxwell-Boltzman DF is «0.2 and at 
which the time-dependent DF is close to a Juttner DF (see Fig|5]), 
corresponds to five Spitzer times or 6.5 x 10 5 yrs for an elec- 
tron number density of 1CT 3 cirT 3 . Using a merging time scale of 
> 10 8 yrs, we found that the value of the critical number density 
< 6.5 x 1CT 6 cirT 3 . This electron number density is much smaller 
than the average gas density in galaxy clusters (see e.g. Sarazin 
1986) and, therefore, it is possible to distinguish a Juttner DF 
from a Maxwell-Boltzman DF in galaxy clusters by the method 
proposed in the Sect. 3. 



5. DFs in galaxy clusters with hard X-ray spectral 
tails 

Galaxy clusters with a hard X-ray excess are promising targets to 
test electron DFs by means of multi-frequency SZ observations. 
This is because these clusters show evidence of either a very 
high temperature plasma or a quasi-thermal emission tail due to 
MHD acceleration mechanisms in the cluster atmosphere. 

Petrosian (2001) estimated the yield in non-thermal 
bremsstrahlung photons and found that a large amount of the 
energy of the non-thermal electrons is transferred to the back- 
ground plasma, so that the ICM should be heated to above its 
observable temperature within 10 Myr. However, it was shown 
by Liang et al. (2002) and Dogiel et al. (2007) that a quasi- 
thermal electron population might overcome this difficulty by 
means of a higher radiative efficiency (and therefore a longer 
overheating time, but see Petrosian & East 2008). Wolfe & 
Melia (2008) also considered a quasi-thermal electron distribu- 
tion when fitting hard X-ray emission, but rather than requiring a 
second-order Fermi acceleration to produce quasi-thermal elec- 
trons, they assumed that quasi-thermal electrons are produced by 
collisions with non-thermal protons. 

A study of the influence of suprathermal electrons on the SZ 
effect was made for the Coma and Abell 2199 clusters by Blasi 
et al. (2000), Shimon & Rephaeli (2002), and Colafrancesco et 
al. (2009). An alternative probe to study the electron distribution 
in galaxy clusters, namely the flux ratio of the emission lines 
caused by FeKa transitions (FeXXV and FeXXVI) was pro- 
posed by Prokhorov et al. (2009). This flux ratio is very sensitive 
to the population of electrons with energies higher than the ion- 
ization potential of a FeXXV ion (which is 8.8 keV). Kaastra et 
al. (2009) demonstrated that the relative intensities of the satel- 
lite lines are sensitive to the presence of suprathermal electrons 
in galaxy clusters. Prokhorov (2009) studies the influence of 
high-energy electron populations on metal abundance estimates 
in galaxy clusters and shows that the effect of high-energy par- 
ticles can be significant. However, the mentioned approaches do 
not allow us to derive electron DFs in galaxy clusters with a hard 
X-ray excess. 

In the section, we propose an approach to derive the electron 
DF and demonstrate that multi-frequency SZ observations are 
promising for this purpose. 

Evidence for non-thermal X-ray emission from the Bullet 
cluster reported by Petrosian et al. (2006) and Million & Allen 
(2009) suggests a possible high-energy subrelativistic electron 
component if these X-ray spectra are interpreted in terms of 
bremsstrahlung emission. Therefore, the Bullet cluster is an in- 
teresting target to test electron DFs. We calculate the SZ inten- 
sity spectra for the Bullet cluster for DFs with and without the 



second electron population proposed by Petrosian et al. (2006). 
These SZ spectra are shown in Figj7] Here we assume that the 
second thermal electron population with k\,T e 2 = 50 keV mim- 
ics the presence of quasi-thermal electrons in the DF and that 
the second electron population is 5% of the first thermal electron 
population. This fraction of high-energy subrelativistic electrons 
is consistent with that derived by Petrosian et al. (2006). For 
the sake of comparison with the results obtained in the previ- 
ous section, we assume that the temperature of the first thermal 
component equals 15.3 keV, which agrees with the temperature 
values of 14. 5^^ keV and 14.8+j^ keV derived by using both 
ASCA and ROSAT (Liang et al. 2000) and by using Chandra 
data (Markevitch et al. 2002), respectively. Note that the total 
DF is normalized. The Comptonization parameter should be in- 
dependently derived from a precise SZ intensity measurement at 
a frequency v = 250 GHz (x=4.4), since the SZ intensities have 
the same value at this frequency for normalized Juttner DFs with 
and without the second electron population (see Fig|7]i. 

For a plasma with these parameters, the decomposition co- 
efficients of the DF with the second electron component approx- 
imately equal to 2.000, 1.207, -0.220, -0.917, -0.741, -0.298 
for the terms 1/2, cos(7r/?), cos(27r/3), cos(3ny3), cos(47iy6), and 
cos(57TjS), respectively. To analyze the DF, we use the method 
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Fig. 7. SZ intensity spectra G(x) = AI(x)/(Ioy) for the cluster 
with the hard X-ray tail for DFs with and without the second 
electron population shown by the dashed and solid lines, respec- 
tively. 

of multi-frequency SZ observations at the frequencies of v=375, 
600, 700, 857 GHz described in the previous section. To study 
our ability to derive DFs from SZ intensity measurements, we 
calculate the differences between the corresponding Fourier co- 
efficients for normalized Juttner DFs with and without the sec- 
ond electron component, which are | AA 1 1 = 0.049, IAA2I 
0.047. |AA 3 | = 0.037, |AA 4 | = 0.045, and |AA 5 | = 0.012. Note 
that the difference | AA5 1 = 0.012 is small and it will be difficult 
to observe this difference by the SZ measurements. However, we 
have checked that this difference has a negligible impact on the 
precision of the approximation of the DF obtained by means of 
Fourier series because of its small values. 
From Eqs.(l23l> and (l25l l we find that the uncertainty in SZ inten- 
sity measurements to distinguish the values of the Fourier coef- 
ficients for normalized Juttner DFs with and without the second 
electron component should be less than: £ = 0.37% to derive the 
A] value, £ = 0.32% to derive the A 2 value, £ = 0.20% to derive 
the A3 value, and £ = 0.1% to derive the A4 value. Therefore, the 
electron DF can be derived by means of the SZ effect if uncer- 
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tainties of SZ observational data are less than 0.1%. 
To show that multi-frequency SZ observations allow us to derive 
the Juttner DF with the second electron component, we calculate 
the SZ intensity (EqO) at the frequencies of y=375, 600, 700, 
857 GHz assuming this DF. We find the decomposition coeffi- 
cients from Eq.(13) by means of the simulated SZ observations 
and show the obtained approximation to the Juttner distribution 
with the second electron population in Figj8]by a dotted line in 
the beta range of ft — [0, 0.7]. The normalized Juttner DFs with 
and without the second electron population are also shown in 
Figj8]by dashed and solid lines, respectively. 




0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 
beta 

Fig. 8. Normalized Juttner DFs with and without the second 
electron population, and the approximation to a Juttner distribu- 
tion with the second electron population obtained by means of 
the simulated SZ observations shown by the solid, dashed, and 
dotted lines, respectively. 

Since the approximate function obtained by means of the 
simulated SZ observations coincides precisely with the normal- 
ized Juttner DF with the second electron population, we con- 
clude that precise observations of the SZ effect can allow to de- 
rive the electron DFs in galaxy clusters with a hard X-ray excess 
if electron DFs are different from a Juttner DF. 

6. SZ intensity maps for the simulated merging 
galaxy cluster 

In the previous sections, we showed that high-frequency spec- 
tral SZ observations at frequencies v = 375 GHz, 600 GHz, 700 
GHz, and 857 GHz provide us with a method to derive DFs of 
electrons in galaxy clusters. 

Here we study the intensity maps of hot clusters at these 
frequencies. Note that all the SZ intensity maps derived in the 
Kompaneets approximation have the same spatial morphology. 
Since the shape of the SZ effect is most sensitive to an electron 
DF at frequencies v = 375 GHz, 600 GHz, 700 GHz, and 857 
GHz (see Sects. 3 and 4), the spatial SZ intensity maps, which 
are derived in the Wright formalism, at these frequencies should 
have different spatial morphology. 

In this section, we show that the spatial morphologies of the 
SZ intensity maps are different at frequencies v = 375 GHz, 600 
GHz, 700 GHz, and 857 GHz. To this aim we use the 3D numer- 
ical hydrodynamic simulations of a merging hot galaxy cluster 
presented in Akahori & Yoshikawa (2010). These authors con- 
sidered a head-on encounter of two free-falling galaxy clusters 
from the turn-around radius (for a review, see Sarazin 2002) . We 
assume in the simulation that the two galaxy clusters have equal 



virial mass of M v ; r = 8 x 10 14 M , and the impact parameter 
is zero (for details on the simulation, see Akahori & Yoshikawa 
2010). 

To produce the SZ intensity maps at frequencies v = 375 
GHz, 600 GHz, 700 GHz, and 857 GHz we use the 3D den- 
sity and temperature maps (see Figs. 2 and 3 from Akahori & 
Yoshikawa 2010) for the simulated merging galaxy cluster at a 
time of t=0.5 Gyr, where t=0 Gyr corresponds to the time of the 
closest approach of the centers of the dark matter halos. We cal- 
culated the SZ intensity maps using the Wright formalism and 
the approach described in Prokhorov et al. (2010a, 2010b). The 
intensity maps of the SZ effect at these frequencies derived from 
the simulated maps of the gas density and temperature are shown 
in Figs. [9j[T0l[Hl andQ~2j respectively. 




x (Mpc) 

Fig. 9. Intensity map I/Iq of the SZ effect at a frequency 375 
GHz derived from the numerical simulation in the framework of 
the Wright formalism. 




x (Mpc) 

Fig. 10. Same as Fig|9]but for a frequency of 600 GHz. 
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Fig. 11. Same as Fig(9]but for a frequency of 700 GHz. 
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Fig. 12. Same as Fig(9]but for a frequency of 857 GHz. 



The morphologies of the SZ intensity simulated maps at v = 
375 GHz, 600 GHz, 700 GHz, and 857 GHz are clearly different. 
This is because the SZ effect from hot galaxy clusters at high 
frequencies is sensitive to the relativistic effects in their electron 
DF and cannot be described in the framework of the Kompaneets 
approximation. Note that the morphologies of the SZ intensity 
simulated maps at frequencies v = 128 GHz and 369 GHz in 
Prokhorov et al. (2010a) are instead more similar because the SZ 
effect from cool galaxy clusters is less sensitive to the relativistic 
effects. 

We find that the SZ intensity simulated map at frequency of 
375 GHz is similar to the SZ intensity simulated map derived in 
the Kompaneets approximation, while the SZ intensity simulated 
map at a frequency of 857 GHz is similar to the SZ intensity sim- 
ulated map at a frequency of 217 GHz where the SZ effect in the 
framework of the Kompaneets approximation is zero. The con- 
tribution of relativistic corrections to the SZ signal at frequency 
of 857 GHz dominates over that derived in the Kompaneets ap- 
proximation. 

We find that the maximum of the SZ intensity increment 
at a frequency of 375 GHz (see Fig. |9]l is at the center of this 
map, while the maximum of the SZ intensity increment at a fre- 



quency of 857 GHz (see Fig. [L21 is in the post-shock regions 
(see the Mach number distribution of ICM in Fig. 2 of Akahori 
& Yoshikawa 2010). This is because the gas temperature in the 
post-shock regions is the highest and SZ relativistic corrections 
are therefore the most significant in these regions. 

Comparing Fig. [9] with Fig. [T2l shows that the SZ intensity 
maps at low and high frequencies look different for a hot merg- 
ing galaxy cluster and, therefore, multi-frequency SZ observa- 
tions with a high spatial resolution are necessary to demonstrate 
the different morphologies of SZ intensity maps at low and high 
frequencies. We conclude that measurements of the SZ intensity 
maps at frequencies v = 375 GHz, 600 GHz, 700 GHz, and 857 
GHz are relevant to study the importance of SZ effect relativistic 
corrections. 

Using the X-ray observations of the Bullet cluster and Abell 
2219 (Million & Allen 2009), which are massive merging galaxy 
clusters, we calculated the SZ intensity maps at low and high 
frequencies and have checked that the morphologies of the SZ 
intensity maps at low and high frequencies are different. We will 
address this problem more specifically in a forthcoming paper. 

7. Observational considerations 

We briefly discuss here the detectability of the electron DF 
through SZ effect observations in the optimal frequency sam- 
ple studied above. 

We find that the difference of the SZ effect signals predicted 
for the Juttner and the Maxwell-Boltzmann DFs at 375, 600, 
700 and 857 GHz can be only marginally appreciated by the 
PLANCK-HFI instrument with the nominal 2-year survey sen- 
sitivity (at 1 cr level). Indeed, while the DF difference in the pre- 
dicted SZ effect intensities can be determined at 375 and 600 
GHz channels at the ~ 3.5<x and ~ 3.4cr level, it is impossible 
to determine it for the highest channels, i.e. 700 and 857 GHz of 
PLANCK-HFI. Note also that the Bullet cluster is almost unre- 
solved for the PLANCK-HFI frequency channels. 
However, the 600 and 857 GHz frequency channels are also cov- 
ered by the HERSCHEL-SPIRE instrument, whose sensitivity 
is sufficient to detect the predicted SZ effect intensity differ- 
ence necessary to distinguish between the Juttner and Maxwell- 
Boltzmann DFs at 600 and 857 GHz. We also notice that the 
Bullet cluster is a fully resolved source for the HERSCHEL- 
SPIRE instrument. 

An important point to be addressed is the impact of the confusion 
noises at these high-frequency channels to distinguish the DF ef- 
fects in the SZ effect observations. Relevant sources of confusion 
for SZ effect observations at high frequencies are the CMB fluc- 
tuation, the unresolved point-like millimeter sources emission, 
and the diffuse emission from the Galaxy. For an experiment 
that can spatially resolve the Bullet cluster, the CMB confusion 
is negligible on sub-arcmin scales, while the unresolved point- 
like source emission and the Galaxy emission increase with fre- 
quency, which provides the major confusion sources. For cluster 
temperatures on the order of 15 keV (as in the Bullet cluster), 
the diffuse Galaxy emission is the dominant confusion noise. 
Therefore, a multi-frequency observational strategy is required 
to properly monitor and subtract the confusion noises that have 
a different frequency behavior with respect to the SZ effect spec- 
trum. Nonetheless, the confusion noise level adds up to the in- 
strument noise to reduce the possibility of the SZ effect obser- 
vations in order to derive the electron DF. Longer exposure with 
spectroscopic instruments operating in this high-frequency range 
are required to reach a precision level on the order of ~ 0. 1 % to 
achieve a good statistical control of the systematics and a good 
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statistical confidence level of the parameters necessary to recon- 
struct the electron DF. 



8. Conclusions 

The SZ effect is an important tool for cosmology and for the 
astrophysical study of clusters of galaxies (for a review, see 
Birkinshaw 1999). It measures the pressure of an electron popu- 
lation integrated along the line of sight as long as free electrons 
are non-relativistic. Relativistic corrections of the SZ effect al- 
low us to measure the plasma temperature. 

Although it has already been noticed that the use of the 
correct relativistic equilibrium distribution is essential for the 
proper interpretation of measurements of the SZ effect, no stud- 
ies of this problem have been performed in detail so far. 
The relativistic kinetic theory, on which the DF derivation is 
based, is still a subject of numerous debates. The relativistic 
analogue of the Maxwell-Boltzman velocity DF has been pro- 
posed by Juttner (1911). However, alternatives to a Juttner DF 
have been discussed by Horwitz et al. (1989) and, recently, by 
Lehmann (2006) and Dunkel & Hanggi (2007). 

In this paper, we showed how to separate the CMB distor- 
tions (caused by the SZ effect) that are caused by a departure 
from the diffusive approximation given by Kompaneets (1957) 
from those that are caused by using a relativistic correct DF in- 
stead of a Maxwell-Boltzman DF. We propose here a method 
based on Fourier analysis to derive a velocity DF of electrons by 
using SZ observations at four frequencies. 

We found that the SZ intensity change owing to using a rel- 
ativistic correct DF instead of a Maxwell-Boltzman DF con- 
tributes a significant part to the total relativistic corrections of 
the SZ effect and is more significant at lower temperatures * 5 
keV than that at higher temperatures a; 15 keV. We conclude that 
the value of the SZ intensity change owing to using a relativis- 
tic correct DF instead of a Maxwell-Boltzman DF will be much 
higher in hot galaxy clusters because the value of the relativistic 

5 /2 

SZ corrections is proportional to T e . 

We proposed a method to derive the DF of electrons us- 
ing SZ multi-frequency observations of massive galaxy clus- 
ters with high plasma temperatures. Using a Fourier analysis 
we found that the approximate electron DF represented by six 
Fourier cosine terms accurately describes the relativistic Juttner 
DF. By means of SZ intensity measurements at four frequencies 
we showed how to derive the approximate DF of electrons. To 
find a suitable sample of four frequencies for deriving the DF of 
electrons, we studied different samples of frequencies. We found 
that the best sample includes high frequencies v=375, 600, 700, 
857 GHz, while the worst sample is in the frequency range 300- 
400 GHz because the matrix Mm, which is determined in Sect. 
3, is ill-conditioned for this frequency range. In the case of the 
frequency sample of v=375, 600, 700, 857 GHz to derive a DF 
of electrons the allowable uncertainties in the SZ intensities are 
about three orders of magnitude larger than those found for the 
frequency sample v = 300, 320, 340, 360 GHz. Using Monte- 
Carlo simulations of SZ observations with the SZ intensity ob- 
servational uncertainty of 0.1%, we showed that it is possible 
to distinguish Juttner and Maxwell-Bolzman DFs by means of 
these SZ observations. Therefore, the SZ effect provides us with 
a promising method to study DFs of electrons in massive galaxy 
clusters that contain hot plasmas with temperatures ^15 keV. 

We considered in our analysis the relaxation of a system 
of electrons with Coulomb interactions and conclude that the 
Juttner distribution function is an appropriate approximation to 



the universal (equilibrium) electron distribution in hot merging 
clusters. 

We applied this method to derive the DF of electrons using 
multi-frequency SZ observations in galaxy clusters where hard 
X-ray tails were detected. We demonstrated the ability of SZ 
multi-frequency observations to derive the electron DF for the 
Bullet cluster and found that a precision of SZ intensity measure- 
ments of »0. 1 % is required. This method is independent of those 
proposed by Blasi et al. (2000), Shimon & Rephaeli (2002), 
and Colafrancesco et al. (2009), and can be used to distinguish 
among different interpretations of the X-ray excess. Although 
our method requires a higher precision of SZ observations com- 
pared with other methods, it permits us to directly study electron 
DFs rather than studying the supplementary electron pressure 
caused by the presence of quasi-thermal particles. This is an im- 
portant advantage of our method. 

Using the 3D hydrodynamic numerical simulations of a hot 
merging galaxy cluster, we also demonstrated that the SZ spa- 
tial intensity maps of the simulated hot merging galaxy cluster 
at frequencies of 375 GHz, 600 GHz, 700 GHz, and 857 GHz 
are different. This is because of the important role of SZ rela- 
tivistic corrections at high frequencies. New SZ multi-frequency 
measurements with a high spatial resolution should confirm our 
conclusion. 

The next generation of SZ effect experimental techniques 
(as those outlined in Colafrancesco & Marchegiani 2010) are 
needed to reach the required high sensitivity with the purpose of 
studying electron DFs by means of multi-frequency SZ observa- 
tions. 
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